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The decay ^c(3900)^ —>■ is searched for using data samples collected with the BESIII detector 

operating at the BEPCII storage ring at center-of-mass energies ^/s = 4.23 and 4.26 GeV. No 
significant signal for the Zc(3900)^ is found, and upper limits at the 90% confidence level on the 
Born cross section for the process e''"e“ —>■ ^c(3900)^7r^ —>■ cu7r"''7r“ are determined to be 0.26 and 
0.18 pb at yfs = 4.23 and 4.26 GeV, respectively. 
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PACS numbers: 14.40.Rt, 13.66.Be, 14.40.Pq, 13.25.Jx 


I. INTRODUCTION 


Recently, in the study of e“'"e“ —> a distinct 

charged structure, named the Zc(3900)^, was observed 
in the spectrum by BESIII [l| and Belle 0. Its 

existence was confirmed shortly thereafter with CLEO-c 
data 3 . The existence of the neutral partner in the decay 
Zc(3900)° —>■ J/'07r° has also been reported in CLEO-c 
data 0 and by BESIII 0. The Zc(3900) is a good can¬ 
didate for an exotic state beyond simple quark models, 
since it contains a cc pair and is also electrically charged. 
Noting that the Zc(3900) has a mass very close to the 
D*D threshold (3875 MeV), BESIII analyzed the pro¬ 
cess —>■ Tr^{DD*)^, and a clear structure in the 

{DD*y^ mass spectrum is seen, called the Zc(3885). The 
measured mass and width are (3883.9±1.5±4.2) MeV/c^ 
and (24.8 ± 3.3 ± II.0) MeV, respectively, and quan¬ 
tum numbers = 1“'' are favored [^. Assuming the 
Zc(3885) —>■ DD* and the Zc{3900) signals 

are from the same source, the ratio of partial widths 
determined to be 6.2 ± 1.1 ± 2.7. 

The observation of the Zc(3900) has stimulated many 
theoretical studies of its nature. Possible interpreta¬ 
tions are tetra-quark Q, hadro-charmonium Q, D*D 
molecule @ and threshold effects (MU. Lattice QCD 
studies provide theoretical support for the existence of 
A(3872) [T^ but not for the ^^(3900) [T3 - [T7l| . However, 
those studies were carried out on small volumes with un- 
physically heavy up and down quarks. It is also worth 
noting that no resonant structure in J/'07r is observed in 
—» J/'07r~*~7r~ by LHCb [l^, in by 

Belle or in 'yp J/'i/^Tr+n by COMPASS [2^. 

The decay properties of a state can provide useful infor¬ 
mation on its internal structure. There are three impor¬ 
tant decay modes for charmonium-like states: (i) “fall- 
apart” decays to open charm mesons; (ii) cascades to hid¬ 
den charm mesons; and (iii) decays to light hadrons via 
intermediate gluons. In addition, as shown in Ref. I,[13: 
an enhancement near the DD* threshold can be pro¬ 
duced via rescattering of hidden or open charm hnal 
states. Decays of the Zc(3900) to light hadrons can play 
a unique role in distinguishing a resonance from thresh¬ 
old effects, because the decay mode with cc annihila¬ 
tion involves neither hidden nor open charm final states. 
However, theory estimates of annihilation widths to light 
hadrons are only order of magnitude due to uncertainties 
of wave function effects and QCD corrections [2l|, 1^ . 
A sizeable Zc(3900) decay width to light hadrons might 
be expected in analogy to rje or XcJ into hadronic hnal 
states. 

Among a large number of hadronic hnal states that are 
available for a /‘^(J^) = 1+(I“'') resonance decay, wtt is 
one of the typical decay modes which are not suppressed 
by any known selection rule. In this paper, we report a 
search for Zc(3900)^ —>■ based on e+e“ annihilation 


samples taken at center-of-mass (CM) energies y/s = 4.23 
and 4.26 GeV. The data samples were collected with the 
BESIII [23 detector operating at the BEPCH storage 
ring. The integrated luminosity of these data samples are 
measured by analyzing the large-angle Bhabha scattering 
events with an uncertainty of 1.0% [l^] and are equal to 
1092 pb“^ and 826 pb“^, for ^/s = 4.23 and 4.26 GeV, 
respectively. 


II. BESIII EXPERIMENT AND MONTE CARLO 
SIMULATION 


The BESIII detector, described in detail in Ref. [^ . 
has a geometrical acceptance of 93% of 47r. A small¬ 
cell helium-based main drift chamber (MDG) provides 
a charged particle momentum resolution of 0.5% at 1 
GeV/c in a 1 T magnetic held, and supplies energy- 
loss {dE/dx) measurements with a resolution of 6% for 
minimum-ionizing pions. The electromagnetic calorime¬ 
ter (EMG) measures photon energies with a resolution of 
2.5% (5%) at 1.0 GeV in the barrel (end-caps). Particle 
identiheation (PID) is provided by a time-of-hight sys¬ 
tem (TOE) with a time resolution of 80 ps (110 ps) for 
the barrel (end-caps). The muon system, located in the 
iron hux return yoke of the magnet, provides 2 cm po¬ 
sition resolution and detects muon tracks with momenta 
greater than 0.5 GeV/c. 

The GEANT4-base d Monte Carlo (MG) simula¬ 
tion software BOOST includes the geometric descrip¬ 
tion of the BESIII detector and a simulation of the 
detector response. It is used to optimize event selec¬ 
tion criteria, estimate backgrounds and evaluate the de¬ 
tection efficiency. We generate signal MC samples of 
e+e” —>■ Zc(3900)^7r^ —cutt+tt” uniformly in phase 
space, where the w decays to 7r+7r“7r°. The decays of 
w —>■ 7r+7r“7r° are gene rated with the OMEGA_DALITZ 
model in evtgen |27l . 1^. I nitial state radiation (ISR) 
is simulated with kkmc p^ . [30l| , where the Born cross 
section of e+e“ — Zc{S900)^tt^ is assumed to follow 
a V(4260) Breit-Wigner (BW) line shape with reso¬ 
nance parameters taken from the Particle Data Group 
(PDG) M- Final state radiation (FSR) effects asso¬ 
ciated with charged particles are handled with PHO¬ 
TOS [23 ■ For studies of possible backgrounds, inclusive 
y(4260) MC samples with luminosity equivalent to the 
experimental data at ^ = 4.23 and ^/s = 4.26 GeV 
are generated, where the main known decay channels are 
generated using evtgen PM with branching frac¬ 
tions taken from the PDG & The remaining events 
associated with charmonium decays are generated with 
LUNDCHARM [33, while continuum hadronic events are 
generated with PYTHIA [13 • QED processes such as 
Bhabha scattering, dimuon and digamma events are gen¬ 
erated with KKMC pi, 113 • 




III. DATA ANALYSIS AND BACKGROUND 
STUDY 
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Tracks of charged particles in BESIII are reconstructed 
from MDC hits. We select tracks with their point of clos¬ 
est approach within ±10 cm of the interaction point in 
the beam direction and within 1 cm in the plane per¬ 
pendicular to the beam. Information from the TOF and 
dE/dx measurements are combined to form PID confi¬ 
dence levels for the tt and K hypotheses; each track is 
assigned to the particle type with the highest confidence 
level. 

Photon candidates are reconstructed by clustering 
EMC crystal energies. The efficiency and energy resolu¬ 
tion are improved by including energy deposits in nearby 
TOF counters. The minimum energy is required to be 
25 MeV for barrel showers (|cos0| < 0.80) and 50 MeV 
for endcap showers (0.86 < |cos0| < 0.92). To exclude 
showers from charged particles, the angle between the 
shower and the extrapolated charged tracks at the EMC 
must be greater than 5°. A requirement on the EMC 
cluster timing with respect to the event start time is ap¬ 
plied to suppress electronic noise and energy deposits un¬ 
related to the event. 

The 7 r° candidates are formed from pairs of photons 
that can be kinematically fitted to the known mass. 
The from this fit with one degree of freedom is required 
to be less than 25. 

Events with exactly four charged tracks identified as 
pions with zero net charge and at least one candi¬ 
date are selected. A five-constraint kinematic fit (5C) is 
performed to the hypothesis of e“'"e“ —>■ Tr+Tr^Tr+Tr^Tr*^ 
(constraints are the 4-moinentuin of the initial e^e“ sys¬ 
tem and the 7 r° mass), and x^c ^ required. If there 
more than one is found in an event, the combination 
with the smallest x^c retained. 

Figure [1] shows the invariant mass distribu¬ 

tion of the TT^TT~TT^ combination with invariant mass clos¬ 
est to the mass of w for the selected candidate e+e“ —^ 
7 r+ 7 r“ 7 r+ 7 r“ 7 r° events at y/s = 4.23 GeV, where promi¬ 
nent 77 , uj and (f) signals are observed. 

uj candidates are selected with the mass window 
|M(7r+^-7rO) closest — Tow I <0.03 GeV/c^, where m,., is 
the nominal mass of the uj taken from the PDG 
Figure [2] shows the M{ujtt^) distribution for the can¬ 
didate events of e+e“ —>■ oitt+tt" at ^/s = 4.23 GeV. 
No sign of a peak near 3.9 GeV/c^ is apparent. The 
shaded histogram in Fig. [2] shows the distribution of 
non-w background for the events in lu sideband regions 
(0.06 < |M(7r+7r“7r°)ciosest - Wwl < 0.09 GeV/c^). 

By studying inclusive MG samples with luminosity 
equivalent to the data at ^/s = 4.23 and 4.26 GeV, the 
background is found to be dominantly from the contin¬ 
uum process e"'"e“ —>■ a; 7 r''' 7 r“. The solid histogram in 
Fig. [5] shows the ujtt^ invariant mass distribution for 
events selected from the inclusive MG sample. 



FIG. 1. The 7r'''7r“7r'’ invariant mass distribution of the 
combination closest to the lo, for the selected e''"e“ 

'K~n~TT^ candidates for the data sample at = 
4.23 GeV. 



FIG. 2. Distribution of M{uj'k^) for the data sample at 
^/s = 4.23 GeV. The dots with error bars are events within the 
u! signal region. The shaded histogram shows events selected 
from the oj sidebands, and the solid histogram shows inclusive 
MG events, which are dominated by continuum events. 

IV. FITTING RESULTS 

We use a one-dimensional, unbinned, extended maxi¬ 
mum likelihood fit to the ujtt^ invariant mass distribution 
to obtain the yield of Zc(3900)^ —>■ ujtt^ events. The sig¬ 
nal probability density function (PDF) is parameterized 
by an iS-wave Breit-Wigner function convolved with a 
Gaussian resolution function and weighted with the de¬ 
tection efficiency: 

where G{M;<t) is a Gaussian function representing the 
mass resolution. The mass resolution of the Zc(3900)^ 
is 1.2 ± 0.1 MeV/c^ at both ^/s = 4.23 and 4.26 GeV, 
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according to MC simulation, p ■ q is the S-wave phase 
space factor, where p is the Zc(3900)^ momentum in 
the e“'"e“ CM frame and q is the w momentum in the 
Zc(3900)^ CM frame. M is the invariant mass of 
and Mq and F are the mass and width of the 2'c(3900)^, 
which are fixed to the results in Ref [l|. e{M) is the 
efficiency curve as a function of the uitt^ invariant mass, 
obtained from signal MC simulation. 

The backg round shape is described by an ARGUS 
function Myjl — (M/mo)'^ ■ exp(c(l — (M/mo)^)), where 
c is left free in the fit and uiq is fixed to the threshold of 

y/s -uIt, [Ill- 

Figure [3lla) shows the fit result for the data sample 
at -y/s = 4.23 GeV. The fit yields 14 ± 11 events for 
the Zc(3900)^ signal. Compared to the fit without the 
Zc(3900)^ signal, the change in InL with A{d.o.f.) = 1 
is 0.74, corresponding to a statistical significance of 1.2a. 
Using the Bayesian method [Slj, Sect.38. 4.1], the upper 
limit for the 2'c(3900)^ signal is set to 33.5 events at the 
90% confidence level (C.L.), where only the statistical 
uncertainty is considered. 

The fit result for the data sample at y/s = 4.26 GeV 
is shown in Fig. [211b). The fit yields 2.2 ± 8.1 events for 
the Zc{3900)^ with a statistical significance of O.lcr. The 
upper limit is 18.8 events at the 90% C.L. 


V. CROSS SECTION UPPER LIMITS AND 
SYSTEMATIC UNCERTAINTY 

The upper limit on the Born cross section at the 90% 
C.L. is calculated as 

cr(e+e“ ^ Zc(3900)=^7rT,Zc(3900)=^ ^ = 

NVL 

_ 1 __ ( 2 ) 

Ant(l + (5)ji3npe(l - ’ 

where is the upper limit on the signal events; £int is 
the integrated luminosity; e is the selection efficiency ob¬ 
tained from signal MC simulation, which are 18.5±0.2% 
and 18.6±0.2% at ^/s = 4.23 and 4.26 GeV, respectively; 
CTe is the systematic uncertainty of the efficiency de¬ 
scribed in next paragraph; is the vacuum polariza¬ 

tion factor obtained by using calculations from Ref. [sl , 
and equal to 1.06 for both energies; (l + <5) is the radiative 
correction factor, equal to 0.844 for ^/s = 4.23 GeV and 
0.848 for = 4.26 GeV obtained using Ref. by 

assuming the line shape of Born cross section a{e'^e~ —>■ 
Zc(3900)^7rT) to be a BW function with the parameters 
of the Y(4260) taken from PDG [H; and B^ and are 
the branc hing fractions of the decay uj —>■ and 

TfO respectively. A conservative estimate of 

the upper limit of the Born cross section is determined 
by lowering the efficiency by one standard deviation of 
the systematic uncertainty. 

The systematic uncertainty of the cross section mea¬ 
surement from Eq. |2| is summarized in Table [H The lu¬ 
minosity is measured using Bhabha events with an un- 




FIG. 3. Results of the unbinned maximum likelihood fit of the 
tJ7r^ mass spectrum of e'^'e” —>■ at (a) = 4.23 GeV 

and (b) ^/s = 4.26 GeV. Dots with error bars are the data. 
The solid curve is the result of the fit described in the text. 
The dotted curve is the .2^(3900)^ signal. The dashed curve 
is the background. 


certainty of 1.0% [ll. The uncertainty in tracking ef¬ 
ficiency for pions is 1.0% per track [|, i.e. 4.0% for the 
track selection in this analysis. The uncertainty in PID 
efficiency for pions is 1.0% per track [|. The uncertainty 
in the photon reconstruction efficiency is less than 1% per 
photon [^ . The uncertainty in the 7r° reconstruction ef¬ 
ficiency is 2.0% [13]. The uncertainty of the kinematic 
fit is estimated by correcting the helix parameters of the 
charged tracks. The detailed procedure to extract the 
correction factors can be found in Ref. [13]. The track 
parameters in MC samples are corrected by these factors, 
and the difference in efficiencies of 0.8% with and with¬ 
out the correction is taken as the systematic uncertainty 
associated with the kinematic fit. An MC sample gener¬ 
ated with Zc(3900)^ —>■ uJTi^ in both S wave and D wave, 
assuming a D/S waves amplitude ratio of 0.1, results in 
a 3% change in detection efficiency. This difference is 
taken as the systematic uncertainty associated with the 
MC production model. The branching ratio value for 
UJ —>■ 7r+7r“7r° comes from the PDG [ll|, and its error is 
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TABLE I. Summary of the relative systematic uncertainties 
of the cross section measurement (in %). 

Source \/s = 4.23 GeV y/s = 4.26 GeV 


Luminosity 

1.0 

1.0 

Tracking 

4.0 

4.0 

PID 

4.0 

4.0 

photon reconstruction 

2.0 

2.0 

7r° reconstruction 

2.0 

2.0 

Kinematic fit 

0.8 

0.8 

Decay model 

3 

3 

Radiative correction 

6 

7 

—>■ 7r^7r“7r°) 

0.8 

0.8 

Total 

9.4 

10.1 


0.8%. In the nominal fit, the radiative correction factor 
and the detection efficiency are determined under the as¬ 
sumption that the production of e+e“ —Zc(3900)^7rT 
follows the y(4260) line shape. Using the line shape 
of CT(e+e“ —Zc(3900)°7r°) measured in Ref. Q as an 
alternative assumption, e(l -I- d) is increased by 6 % for 
^/s = 4.23 GeV and 7% for y/s = 4.26 GeV . The change 
in e(l -h 5) is taken as a systematic uncertainty. The un- 
certainty of the vacuum polarization factor is taken from 
Ref. [^, and is negligible compared with other uncer¬ 
tainties. Assuming that all sources of systematic uncer¬ 
tainties are independent, the total errors are given by the 
quadratic sums. 

To estimate the systematic uncertainties due to the fit 
procedure, we fit under different scenarios, and the upper 
limits obtained at the 90% C.L. for the Z'c(3900)^ sig¬ 
nal yield are summarized in Table HU The effect on the 
signal yield from the fit range is obtained by varying the 
fit range by ±0.1 GeV/c^. The effect due to the choice 
of the background shape is estimated by changing the 
background shape from the ARGUS function to a second 
order polynomial (where the parameters of the polyno¬ 
mial are allowed to vary and the fit range is limited to 
[3.4, 4.08] GeV/c^). The effect due to the resonance pa¬ 
rameters of the 2'c(3900)^ is estimated by varying the 
resonance parameters according to the results in Ref 
The effect due to the mass resolution is estimated by 
increasing the resolution by 10 % according to the com¬ 
parison between the data and MG. The effect due to the 
mass-dependent efficiency curve is estimated by chang¬ 
ing the efficiency curve to a constant function. We take 
the largest number of .Zc(3900)^ events in the different 
scenarios as a conservative estimate of the upper limit: 
V 4230 = 38.0, V 4260 = 18.8. The resulting upper limits 
of the Born cross sections at yfs = 4.23 and 4.26 GeV 
are determined to be 0.26 and 0.18 pb at the 90% G.L., 
respectively. 


VI. SUMMARY AND DISCUSSION 

In summary, based on data samples of 1092 pb“^ 
at yfs = 4.23 GeV and 826 pb“^ at yfs = 4.26 GeV 


TABLE 11. Results of upper limits on the Zc(3900) signal 
yield with various fit procedures. 


Source 

Vi = 4.23 GeV Vi 

= 4.26 GeV 

Fit range 

31.5 

18.5 

Background shape 

38.0 

16.1 

Zc{3900) mass and width 

22.6 

12.2 

Mass resolution 

33.5 

18.8 

Efficiency curve 

33.3 

18.8 


collected with the BESIII detector operating at the 
BEPGII storage ring, a search is performed for the decay 
Zc{3900)^ —>■ WTT^ in e“''e“ —>■ wtt+tt”. No Zc(3900)^ 
signal is observed. The corresponding upper limits on 
the Born cross section are set to be 0.26 and 0.18 pb 
at y/s = 4.23 and 4.26 GeV, respectively. If we assume 
that the Zc(3900)^ observed in e+e“ —>■ J/' 07 r+ 7 r“ [I[ 
and Zc(3885)^ in e+e“ — {DD*)^Tr^ @ are the same 
particle, the decay width of Zc(3900)^ —>• a; 7 r=*= is esti¬ 
mated to be smaller than 0.2% of the ^^(3900)^ total 
width. As LOTT is a typical light hadron decay mode of 
a = 1 +( 1 +) resonance, the non-observation of 

Zc(3900)^ —?► may indicate that the annihilation of 

cc in Zc{3900)^ is suppressed. Gomplementary to the 
searches for Zc(3900) production [isl - l^ . exploring new 
Zc(3900) decay modes may provide a significant input to 
clarify its dynamical origin. 
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